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Abstract 

We present a comprehensive study of pair production and decay of stops, sbottoms, 
and staus in e + e~ annihilation at LEP2. We give numerical predictions within the Mini- 
mal Supersymmetric Standard Model for cross sections and decay rates, and discuss the 
important signatures. In the case of stau production we also study the polarization of the 
r in the decays fi — ► t\\ 2- 



1 Introduction 



Super symmetry (SUSY) |J requires the existence of two scalar partners and q~R (squarks) 
for every quark corresponding to its two helicity states. For the same reason every lepton I has 
as super symmetric partners the sleptons £l and £r. Whereas in the case of the first and second 
generation the q^, q^ and £l, £r states are to a good approximation also the mass eigenstates, 
this is not expected for the third generation. Quite generally, q~L and qn (£l and £r) mix, the 
size of the mixing being proportional to the mass of the quark q (lepton £). Thus it may be that 
one of the mass eigenstates of the stops, ti, is the lightest squark. It could even be possible that 
it is the lightest visible SUSY particle. If the parameter tan /3, which also enters the mixing, is 
large (tan/? > 10), also the sbottom b± or the stau fi can be relatively light. Therefore, it is 
interesting to study stop, sbottom, and stau production in the energy range of LEP2, that is 
mz < \/s < 192.5 GeV. The present experimental bound from LEP1 for the masses of charged 
super symmetric particles is m > 45 GeV [fj, |3[]. Stronger limits for the stop mass, up to 55 GeV, 
are reported from LEP at \/s = 130 — 140 GeV ||. The D0 experiment at FNAL obtained 
additional mass bounds for the stop |5|, [| excluding the mass range 40 GeV < < 100 
GeV if the mass difference (m^ — m^o) > 30 GeV, where the neutralino Xi is the lightest 
super symmetric particle (LSP). 

Here we shall calculate the production rates for e + e~ — ► t\ti, e + e~ — > b\b\, and e + e _ — > 
f\ f i in the whole accessible mass range of LEP2 as a function of the mixing angle. SUSY-QCD 
corrections as well as initial state radiation (ISR) will be included. We also discuss in detail the 
decay patterns of these particles. The phenomenology of the decays of t, b, and f is different 
from that of the other squarks and sleptons due to their non-negligible Yukawa couplings. 

The framework of our calculations is the Minimal Supersymmetric Standard Model (MSSM) 
0. The parameters which determine the phenomenology of stops, sbottoms, and staus are 
M, M', the soft-breaking SU(2) and U(l) gaugino masses, /x, the higgsino mass parameter, 
tan j3 = v 2/Vl, where v i and v 2 are the vacuum expectation values of the two Higgs doublets, m^, 
the mass of the pseudoscalar Higgs boson A , and Mp (F = Q, U, D, L, E) and Af (/ = t, b, r), 
the soft-breaking parameters which enter the mass matrices of the stops, sbottoms, and staus. 
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We shall use the GUT relations M' = f M tan 2 9 W ~ 0.5M and m$ = — iW sin 2 # w ~ 0.3M, 
where is the gluino mass. 

In the energy range of LEP2 the most important decay modes of £1, b\, and f\ are £1 — > 
ex?, fcXi^, &i — > &X2> an d ^1 ~~ *■ r X?25 ^tXi (assuming m fl < m~i). Here is the lighter 
of the two charginos present in the MSSM and Xi is t ne lightest of the four neutralinos (with 
m^o < m^o < m^o < Tn^o). It is important to note that the masses and couplings of the 
charginos and neutralinos only depend on M, fi, and tan (3. 

In the next section we present the formulae for the sfermion mixing, the production cross 
sections and the decays of £1, bi, and fx, taking into account the mixing and the Yukawa 
couplings. In section 3 we give numerical results for the cross sections of e + e~ — > t\ ti, b\ bi, fx fx 
for various masses as a function of the mixing angle, a numerical analysis of the decays of ti, bi, 
and fx into charginos and neutralinos, as well as a discussion of the main signatures. Section 4 
contains a short summary and conclusions. 

2 Production cross section and decay formulae 

The mass matrix for stops, sbottoms and staus in the (j^, basis (with / — t, b, f) has the 
following form ||: 

M)=( m ^ a T f ) (1) 
/ ^ a f m f mj R J 

with 

m] L = M| + mf + m| cos 2(3 (| - § sin 2 9 W ), mf R = M? + m 2 t - |m| cos 2(3 sin 2 9 W , 
m\ L = M^+m 2 b -m 2 z cos2(3(\-\sin 2 9 w ), m 2 R = + m 2 + |m| cos 2(3 sin 2 9 W , (2) 
m\ h = M| + m 2 - m| cos 2/5 (| - sin 2 m\ R = M| + m 2 + m| cos 2(3 sin 2 

and 

a t m 4 = m t (A t — /i cot (3), a b m b = m b (A b — fi tan (3) , a r m T = m T (A T — \i tan (3) 

The mass eigenstates fi and f 2 are related to fi and by: 

(h\ = ( cos 9 f sin 9 f \ ( f L \ 

\f 2 ) { -sin 9 f cos9 f ) \f R ) ™ 
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with the eigenvalues 



m h = 2 H + m L) =f 2\/ (m /,- m L )2+4 



(5) 



The mixing angle 9i is given by 



cos 9i 



l f m f 



(m 2 ? — w?~ ) 2 + a 2 f m 2 f ' 
v jl h 11 



sin 9, 



(m 2 - — in 2 , ) 2 



\ (m 2 ^ — m 2 ^) 2 + ajm 2 ' 
Hence, | cos Of] > l/v2 if mj L < m? and | cos6^-| < l/y/2 if < . 



(6) 



The reaction e + e — > f\ / 1 proceeds via 7 and Z° exchange. The tree level cross section at 
a center-of-mass energy of yfs is given by: 



(7 



tree 



na 2 N c 
3s 



16g 4 c 4 * ~ 2g2 c / s{s-m z) 



1 



s — m|) 2 + T 2 z m 2 z 



(7) 



where = 1 — c 2 ^ = sin 2 9w, v e = 2 sin 2 9w — 2 anc ^ a e = — §• ^c* is a colour factor which is 
3 for squarks and 1 for sleptons. The Z° coupling to j\ j\ is proportional v j i = 2 (7| cos 2 9f — 
Qfsm 2 9w)- Here Ij and Qf are the third component of the weak isospin and the charge of 
the fermion / (Q e = —1). <j tree shows the typical /3 3 suppression where (3 = (l — 4m%/s) 
is the velocity of the outgoing scalar particles. The interference of the 7 and Z° exchange 
contributions leads to a characteristic minimum of the cross section at 



cos 2 # 



/ 1 min 



Ql 

T 3 



sin 2 9 W 



1 1 /I S \ 2 n -^e + Re 

1 + (1 ~^ )cos 6w WTr1 



where L e = sin 2 Ow — \ and R e = sin 2 ^- The angular distribution has the familiar sin 2, # 



shape, with $ the scattering angle: 



da 



tree 



d cos-# 



sin da 



tree 



(9) 



In the case of squarks QCD radiative corrections are important. The conventional QCD correc- 
tions were calculated in 0, || including the radiation of soft and hard gluons (in 0(a s )). The 
QCD corrections within the MSSM including virtual gluino and squark exchange were com- 
puted in |J (The corrections due to the four-squark interaction is zero in the renormalization 



scheme used). In our numerical calculations we have included both the gluonic corrections and 
those due to gluino and squark exchange. The corrections due to the exchange of supersym- 
metric particles are, however, small in the energy range of LEP2. Moreover, we have taken into 
account initial state radiation (ISR) (n 



The sfermion interaction with neutralinos and charginos is given by ]TT 



C = gf (a{ k P L + b{ k P R ) xlfi + g f (Z£ P L + fc£ P R ) x- U + hx. (10) 



with 



It 


= —Vji cos 6i + Y t Vj 2 sin 9+, 


ut 


= Y b U j2 cos 0- t , 


Jt 
'2.7 


= Vji sin 6i + Y t Vj 2 cos 9$, 


ut 

K 2 j 


= - Yb U j2 sin h 


l b 

\i 


= -Uji cos 9 b + Y b U j2 sin 9 b , 


Lb 
lj 


= Y t V j2 cos b , 


'2.7 


= Uji sin 9- b + Y b U j2 cos 9~ b , 


k b 


= -Y t V j2 smO b , 


]T 

''1.7 


= —Uji cos Of + Y T Uj 2 sin Of, 


h T 
K lj 


= 0, 


JT 


= Uji sin 0~ b + Y T Uj2 cos Of, 


K 2 j 


= 0, 



(11) 

(12) 
(13) 



( a {k\ = ( cos Of sin f \ fh[ k \ fb{ k \ = f cos0 f smO f \ f f[ k \ 

h\ k = Y t (N k3 sin - N M cos 0) , j\ k = sin w iV fel - ^ ( i - § sin 2 9 W ) £^ , 

h Rk = Y t (N k3 sin (3 - N k4 cos (3), f Rk = -^ sin 9 W (tan %7V fe2 - JV fc i) , 

= ~ Y b (Nk3 cos /3 + A'/, | sin /?) , = ^ sin %iV fea + (| - ± sin 2 %) 
= ~ Y b ( N k3 cos /? + N k4 sin /5) , f b Rk = ^ sin Mtan %^2 - N kl ), 

h T Lk = -Y T (N k3 cos /3 + N k4 sin (3), f£ k = V2 sin 9 w N kx + ^2 (± - sin 2 6w)£fc, 
h T Rk = -Y T (N k3 cos /3 + N k4 sixi (3), f T Rk = V2 sin 9 W (tan 9 w N k2 - N kl ), 

for the sfermion-fermion-neutralino interaction, and 

(15) 

(16) 
(17) 

for the sfermion-fermion-chargino interaction. Nij is the 4x4 unitary matrix diagonalizing the 
neutral gaugino-higgsino mass matrix in the basis 7, Z , H® cos (3 — sin (3, H® sin (3 + cos /3 
T2\. Uij and Vij are the 2x2 unitary matrices diagonalizing the charged gaugino-higgsino mass 



matrix |13|]. We choose a phase convention in which N^, U^, and Vy are real. Yf denotes the 



Yukawa coupling, 

Y t = m t /(V2mw sin/?), Y b = m b /(V2 m w cos/5), Y T = m T / (V2 m w cos (3) . (18) 
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The respective decay widths then are 



r(/< - fxl) 



g 2 \z(m 2 f , m 2 f ,m? ) 

Ji J X-k 



167rm| 

Ji 



(a 2 ik + 6i fc )(m| -m 2 f - m|o) - Aa ik b ik m f m^ (19) 



and 



g 2 \ 1 2(m 2 r , m 2 fl ,m 2 ±) 



Ji 



(l 2 j + kfyim 2 -. - m 2 , - ml±) - ^l i jk ij m r m ~± 



where X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. 



(20) 



3 Numerical results 

3.1 Stop ti 

The total cross sections for the process e + e~ — > t\t\ at y/s = 175 GeV and y/s = 192.5 GeV 
are shown in Fig. | as a function of | cos9^\ for several mass values of t%. At y/s = 192.5 
(175) GeV, for a stop mass of 80 GeV the cross section reaches 0.35 (0.18) pb. Assuming an 
integrated luminosity of 300 (500) pb^at yjl = 192.5 (175) GeV one can thus expect ~ 60 
to 105 (50 to 90) t\ti events for ~ 80 GeV. Moreover, the cross section shows a clear 
dependence on the mixing angle for mf < 80 GeV and |cos^| > 0.6. In this region cross 
section measurements should therefore allow to determine the stop mixing angle once the mass 
of ti is known. The importance of radiative corrections is illustrated in Fig. where we show 
the conventional QCD corrections in 0(a s ), the corrections due to gluino exchange as well as 
the ISR corrections at = 192.5 GeV as a function of mi for cos6^ = 0.7. The gluonic 
corrections enhance the tree level cross section rising from 17 to 41% for stop masses in the 
range of 45 to 85 GeV. The corrections due to the exchange of supersymmetric particles are 
< 1.2% for rrig = 200 GeV and m 4 ~ 2 = 250 GeV and depend on the stop mixing angle. Initial 
state radiation turns out to alter the tree level cross section from ~ 1.2% to ~ —21% for 
= 45 to 85 GeV. 

Assuming m~ tx < mi~ the main decay modes of t\ are t\ — > CX\ and t\ — > bxt- The latter 
decay has practically 100% branching ratio if it is kinematically allowed. As xt further decays 
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into X? l4 or Xi the signature is two acoplanar b jets accompanied by two leptons + large 
missing energy (1$), or single lepton + jets + 1$, or jets + If. Here the b tagging technique can 
be used to extract the signal. Moreover, in this case the xt wm m ost likely be observed first 
and the information from its decay properties can then be used to identify the t\. On the other 
hand, if < m~+ + m& the flavour changing decay t\ ^ cxj has practically 100% branching 
ratio. The signature is then two acoplanar jets + 1$. Generally, for t\ — > cxi the invisible 
energy is larger than in case of t\ — > bxt- In Fig. |3] we show the domains of the ti decay 
modes in the (M, (/,) plane for = 80 GeV and tan (3 = 2 (there is a small stripe where the 
decay t\ — > cx\ is a l so possible). If, however, m~± + m h > > + m b (+m e ) the decays 
t\ — > bvet or t\ — > b£u£, proceeding via a virtual xt can compete with the decay into cx±- In 
this case the signature is again 2b + 2£ + If. 

If the lifetime of t\ is longer than the typical hadronization time of (9(10~ 23 s), i. e. T < 0.2 
GeV, t\ hadronizes first into a colourless (tig) or (tiqq) bound state before decaying. This 
is generally expected in case of t\ — > CX\ and t\ — > bugi, blue since these decays involve the 



electroweak coupling twice ||14|| . However, also the width of the t\ decay into bxt can be smaller 
than the hadronization scale as illustrated in Fig. f|. Here we show the width of ti — > bxt as 
a function of cos^ for m% x = 85 GeV, m-± ~ 60 GeV, tan/3 = 2, and three scenarios with 
M -C |/i| <C M, and M ~ \fi\ respectively. 

If, as for the curves (a) and (b), M is much smaller than the lighter chargino is gaugino-like. 
In this case the ti-b-Xi interaction is dominated by Vn cosdi in eq. (p!5|) for |cos^| > 0.3. 
The deviation of the decay width from the cos 2 Q~ t shape and the asymmetry in the sign of \x 
are due to constructive and destructive interferences with the term proportional to the top 
Yukawa coupling (Y t Vi 2 sin 0i) which becomes especially important for t\ ~ t R and increases 
with decreasing 

On the other hand, a higgsino-like xt strongly couples to £r, the dominant coupling being 
proportional to the top Yukawa coupling (|/i| -C M, curves (c) and (d)). Thus the decay width 
is large if | cos# t ~| ~ but goes to zero for t\ ~ tz- 

For M ~ which is shown in (e) and (f), a complicated interplay of gaugino and higgsino 
couplings gives rise to an intricate dependence on the mixing angle and a large asymmetry in 
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the sign of \i. 

The process of fragmentation of t\ was discussed in detail in ref. f|. Fast moving stops 
first radiate off gluons at small angles. This process can be treated perturbatively. After that 
the (non-perturbative) hadronization phase follows by forming (tiq) and (t^qq) hadrons. For 
/3 ~ 1/2 the energy loss of t\ due to gluon radiation and due to hadronization is of comparable 
size. Near the threshold, the gluon emission suffers a /3 4 suppression. 

Monte Carlo studies for t\ pair production have been performed within the CERN-LEP2 
Workshop 1995 [fT5|| . They have mainly concentrated on the decay t% ^ ex? since in this case 
t\ would most probably be the first SUSY particle to be discovered. For the simulation of t\ 
fragmentation different approaches have been used. The conclusions have been a 5<r discovery 
reach of m ?1 ~ 75 to 90 GeV and a 95% CL exclusion of m ix ~ 84 to 92 GeV at yfs = 190 GeV, 
for C = 300 pb _1 depending on the stop mixing angle and the mass of the lightest neutralino. 
In case of t% — > b\t the experimental reach for is ~ 85 GeV. 

3.2 Sbottom h 

According to eq. (||) a considerable b^-bji mixing is possible if tan/? is large (tan/3 > 10). In 
this case the sbottom b\ can be rather light. Hence, it is also interesting to discuss b\ search at 
LEP2. 

The total (SUSY-QCD & ISR corrected) cross sections of bi pair production at ^Js = 175 
GeV and = 192.5 GeV are shown in Fig. |5] as a function of | cos^l for several mass values 
of b\. As can be seen, the dependence on the mixing angle is even more pronounced than in 
case of ti production. The cross section of e + e~ — > b\ b\ is smaller than the t\ t\ cross section 
by a factor of 4 to ~ 5/4 for |cos^| in the range of to 1. Again, radiative corrections 
are important. Standard QCD and ISR corrections are similar to the stop case, i. e. gluon 
corrections enhancing the tree level cross section from 17 to 41% for b\ masses in the range 
of 45 to 85 GeV and initial state radiation lowering it by up to 21%. The corrections due to 
gluino exchange are < 1% for m~ g = 200 GeV and m~ b2 = 250 GeV. 
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Assuming m bi < m§ the main decay modes of b\ are b\ — > bxi and b± — > b\2i the second 
decay being possible in the parameter region approximately given by M < m b — m b or |/x| < 
77tg — m&. For the 6%? channel the signature is two acoplanar b jets + missing energy 1$. If the 
b\ decays into 6%° the b jets are acompanied by additional jets and/or leptons from x° — * Xi QQ 
and/or x\ —>■ Xi^- b tagging will help to enhance the signal. The domains of the b\ decays in 
the (M, /x) plane are shown in Fig. |6] for m bi = 80 GeV and tan (3 = 40. 

In Fig. [7] we show the branching ratio for b\ — ► b Xi as a function of cos 9 b for m bi = 80 GeV, 
tan/? = 40, |/x| = 500 GeV, and M = 55, 65, and 75 GeV. As can be seen, the branching ratio 
is highly dependent on the mixing angle, if M <C |/x| and both Xi and X2 are light. In this case 
bi mainly decays into bxi if | cos^^| ~ 0, whereas for b\ ~ 6^ the decay into 6x2 dominates. 
Evidently, the dependence on the mixing angle weakens with increasing x® masses. If, however, 
M ~ or M ^> the decay b\ — > 6 i s the dominant mode. The /x dependence of the 
branching ratio is shown in Fig. H for = 80 GeV, | cos^| = 0.72, tan/3 = 40, and several 
values of M. As the parameter Af, should not become unnaturally large in supergravity models 
jl6f , for large values of tan/3, the parameter /x determines the sign of cos^ (see eq. (^|) and 
eq. (|D). For this reason, we have chosen cos 9 b and /x such that they have the same sign in 
Fig. 0, Fig. g and Fig. |. 

As in the t\ case the decay width of b\ can be significantly smaller than 0.2 GeV as illustrated 
in Fig. p. Here we show the total width of 61 as a function of /x for m~ bi = 80 GeV, | cos 9 b \ = 0.72, 
tan/3 = 40, and several values of M. The decay width only exceeds the hadronization scale 
(r > 0.2 GeV) if Xi is very light ( < 30 GeV) or if the light neutralinos have strong higgsino 
components (small |/x|). It clearly increases with decreasing |/x| as then X12 become lighter 
and the higgsino couplings proportional to Yj, gain importance. So does the dependence on the 
mixing angle, which shows minima at cos^ = ±1 and cos 9 b = where one decay channel (&Xi 
or 6x2) dominates. 



DELPHI has studied b\ search when b\ decays into bxi [ 15 [ . Their conclusion was that in 
this case the discovery potential for b\ is similar to the one of t% — > cxi, i- e. m b < 75 to 90 
GeV depending on 9 b and m^o. 
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3.3 Stau fi 



As in the sbottom sector one expects large t^-tr mixing for high values of tan (3. This would 
lead to a light f\ which could also lie in the energy range of LEP2. 



The f\ T\ production cross sections are plotted in Fig. |10] at yfs = 175 GeV and y/s = 192.5 
GeV as a function of | cos#f | for several f\ masses. As can be seen, the dependence on the 
mixing angle is much weaker for f\ f\ than for t\ t\ and b\ b\ production. At yfs = 192.5 GeV 
the cross section is of 0(0.1) pb for = 80 GeV and goes up to 0.55 pb for a fi mass of 50 
GeV. This corresponds to a production rate of ~ 30 (25) events at y/s = 192.5 (175) GeV for 
an integrated luminosity of C = 300 (500) pb -1 and m fl = 80 GeV. Of course, here only ISR 
corretions have to be taken into account. 

The simplest signature of fi f\ production is 2r+^ with the r leptons coming from f\ — > r Xi 
decays. If the decay modes into r x° an d v T Xi are al so kinematically accessible the signature 
is 2r + jets and/or leptons + 1$, or single r + jets and/or leptons + 1$, or jets and/or leptons 
+ 1$ due to cascade decays of x° and Xi- The latter case occurs when both staus decay via 
fx — > v T Xi (fi — > v T xt)- The parameter domains of the various f\ decays in the (M, \i) plane 
are shown in Fig. |11] for = 80 GeV and tan/? = 40. If e, ft, and/or v are lighter than f\ 
then also three-body decays into these particles plus neutrino(s), lepton(s), or quark pairs are 
possible. 



In Fig. |12| we show the branching ratios of f\ decays into r Xi, t an d v T Xi as a function 
of cosflf for m fl = 80 GeV, tan/3 = 40, (a) M < and (b) < M. As in section 3.2, 
we have chosen cos Of and fi such that they have the same sign in order to avoid unnaturally 
large values of A T . In case (a), i. e. M <C |//|, if cos^ < 0.5 the decay f\ ^ rx? has > 80% 
branching ratio as tr does not couple to a gaugino-like x% or Xi ■ O n the other hand, the decays 
into r x® an d v T Xi play an important role for f\ ~ Tl, i- e. cos^^ ~ 1. In contrast to that, for 
higgsino-like light neutralinos and charginos, i. e. <Mas shown in (b), BR(fi — > txD 
and BR(fi — > v T \\) are °f comparable size for f\ ~ tr, whereas for fi ~ ?l the fi-v T -Xi 
coupling vanishes and the decay into TXi has ~ 100% branching ratio. The fi — > decay 
mode is negligible in this case. 
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For staus a new interesting aspect comes into play: Due to the supersymmetric versions 
of gauge and Yukawa interactions the sfermion-fermion-gaugino interaction conserves chirality 
while the sfermion-fermion-higgsino interaction flips it. The polarization of the final state 
fermion thus depends on the sfermion (L-R) mixing as well as on the gaugino-higgsino mixing. 
As, in contrast to other leptons, taus decay in the detector, one can determine their average 
polarization through the energy distribution of their decay products [|17]]. Hence, one has an 
additional opportunity to obtain information on the t^-tr and the gaugino-higgsino mixing by 
measuring the average polarization of tau leptons coming from f\ — > T\° k decays [|18|j . 
The average polarization of the r leptons is given by: 

~ ( x = Ek BR(f! -> r R xl) - BR(f! -> r L xl) ( , 

[T) Efc BR(fi — > Tr xl) + BR(f 1 — > t l xl) { ] 

V(t) depends on the stau mixing angle Of and on the mixing matrix Ny of the neutral gaugino- 
higgsino sector. If the lighter neutralinos are gaugino-like only Nki and N^ 2 give sizeable 
contributions. Thus, in this case V(t) ~ +1 for | cos#f | ~ 0, whereas V[f) ~ —1 for | cos#f | ~ 
1. On the other hand, in case of higgsino-like x\ and one finds Vir) ~ —1 for | cos^f | ~ 
and V(t) ~ +1 for |cos^f| ~ 1. This is illustrated in Fig. 12 where we show the average 
polarization of tau leptons arising from f\ — > r x\ decays as a function of cos Of for = 80 
GeV, tan/3 = 40, and three scenarios of (a) M -C (b) M, and (c) M ~ In case 

of M ~ where the lighter neutralinos have both, gaugino and higgsino, properties V(r) ~ \ 
and varies only little with the stau mixing angle. 



OPAL has studied the experimental aspects of a stau search at LEP2 for mf R [ ^5| . 
At y/s = 190 GeV and for £ = 300 pb" 1 they obtained a 5cr detectability for stau pairs of 
mf R ~ 70 to 83 GeV for neutralino masses in the range of 20 to 72 GeV. However, they have 
neglected the interesting possibility of f L -f R mixing. 

4 Conclusions 

We have discussed the phenomenology of stop, sbottom, and stau pair production and decays 
at LEP2 paying particular attention to the sfermion L-R mixing. Analytical formulae have 



10 



been given for the sfermion mixing, for the cross sections of t\, bi, and fi pair production 
in e + e~ annihilation, and for the widths of the decays of these particles into neutralinos and 
charginos. We have presented numerical predictions for these production and decay processes 
and analyzed their dependence on the SUSY parameters. It has turned out that due to L-R 
mixing and non-negligible Yukawa couplings the phenomenology of stops, sbottoms, and staus 
can be significantly different from that of first and second generation sfermions. Moreover, we 
have shown that in case of t\ and 61 hadronization effects can be important. For f x we have 
also discussed the dependence of the average polarization of taus arising from f\ — > r decays 
on r L -f R and gaugino-higgsino mixing. 
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Figure captions 



Figure 1: Total cross section for e + e — > fi in pb at y^i = 175 GeV (dashed lines) and 
y/s = 192.5 GeV (solid lines) as a function of the stop mixing angle for t\ masses of 50, 60, 70, 
80, and 90 GeV. 

Figure 2: Gluon, gluino-stop, initial state, and total radiative corrections relative to the tree 
level cross secion of e + e~ — > t\ t\ at = 192.5 GeV as a function of for cosOf = 0.7. 

Figure 3: Parameter domains in the (M, /x) plane for the various i\ decay modes, for m~ ti = 
80 GeV and tan/5 = 2. The grey area is excluded by LEP1. 

Figure 4: Decay width of i x — > &x+ for = 85 GeV, m-- ~ 60 GeV, tan/3 = 2, (a) M — 51.2 
GeV, /i = -500 GeV, (b) M = 72.0 GeV, /i = +500 GeV, 1 (c) M = 400 GeV, n = -49.4 GeV, 
(d) M = 400 GeV, n = +75.5 GeV, (e) M = 50 GeV, /j = -47.9 GeV, and (f) M = 135 GeV, 
li = +132 GeV. 

Figure 5: Total cross section for e + e~ — > foi &i in pb at — 175 GeV (dashed lines) and 
yfs = 192.5 GeV (solid lines) as a function of the sbottom mixing angle for b± masses of 50, 60, 
70, 80, and 90 GeV. 

Figure 6: Parameter domains in the (M, /i) plane for the b\ decay modes for m~ bi = 80 GeV 
and tan/5 = 40, (a) b\ — > bxl and (b) b\ — > 6x? 2 - The grey area is excluded by LEP1. 

Figure 7: Branching ratio of b\ — > 6x? i n percent as a function of the sbottom mixing angle 
for mi; = 80 GeV, tan/3 = 40, = 500 GeV, and M = 55, 65, and 75 GeV. 



Figure 8: Branching ratio of &i — > 6x? i n percent as a function of [i for = 80 GeV, 
cos0j = 0.72, tan/3 = 40, and M = 55, 65, and 75 GeV. 



Figure 9: Total decay width of b\ in GeV as a function of [i for m~ bi = 80 GeV, cosflg = 0.72, 
tan/3 = 40, and M = 55, 65, 75, and 150 GeV. 

Figure 10: Total cross section for e + e~ — > fifi in pb at y/s = 175 GeV (dashed lines) and 
\fs = 192.5 GeV (solid lines) as a function of the stau mixing angle for fx masses of 50, 60, 70, 
80, and 90 GeV. 
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Figure 11: Parameter domains in the (M, /i) plane for the various t\ decay modes, for m fl = 
80 GeV and tan/3 = 40 with (a) n -> r*?, (b) fi -> r^S, , and (c) fi -> r^?, "tXi 
The grey area is excluded by LEP1. 

Figure 12: Branching ratio of f± — > TXi (solid lines), T\ — > rf" (dashed lines), and fx — > z/ T xJf 
(dashdotted lines) in percent as a function of the stau mixing angle for m fl = 80 GeV, tan /3 = 
40, (a) M = 65 GeV, = 500 GeV, and (b) M = 300 GeV, = 60 GeV. 

Figure 13: Average polarization of r leptons arising from f\ — > decays as a function of 
the stau mixing angle for m fl = 80 GeV, tan/9 = 40, (a) M = 65 GeV, = 500 GeV (solid 
line), (b) M = 300 GeV, = 60 GeV (dashdotted line), and (c) M = 100 GeV, = 100 
GeV (dashed line). 
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